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SUMMARY 


An analysis has been made of atmosphere entries for which the 
vehicle lift -drag ratio was modulated to maintain specified maximum 
decelerations and/or maximum deceleration rates. The part of the 
vehicle drag polar used during modulation was from maximum lift coeffi- 
cient to minimum drag coefficient. The entries were at parabolic velocity 
and the vehicle maximum lift -drag ratio was 0.5. Two-dimensional trajec- 
tory calculations were made for a nonrotating, spherical earth with an 
exponential atmosphere. The results of the analysis indicate that for a 
given initial flight -path angle, modulation generally resulted in a reduc- 
tion of the maximum deceleration to 60 percent of the unmodulated value or 
a reduction of maximum deceleration rate to less than 50 percent of the 
unmodulated rate. These results were equivalent, for a maximum decelera- 
tion of lOg, to lowering the undershoot boundary 2k miles with a res ult ing 
decrease in total convective heating to the stagnation point of 22 percent. 
However, the maximum convective heating rate was increased l8 percent; the 
maximum radiative heating rate and total radiative heating were each 
increased about 10 percent . 


INTRODUCTION 


Man’s performance of useful duties during atmosphere entry requires 
that the decelerations and the deceleration rates be kept within prescribed 
limits. For nonlifting vehicles, the limiting of deceleration and deceler- 
ation rates by modulation of drag has been treated in references 1, 2, 
and 3- For lifting vehicles, limited studies of the effects of mod ul ation 
of lift and drag to limit the deceleration were made in references 4 and 5. 
More extensive studies were made in references 6 and 7 but, a dmi ttedly, 
the reductions in maximum deceleration afforded by the use of modulation 
were conservative, especially for vehicles with maximum lift-drag ratio 
as low as 0.5, because of the assumptions on which the analyses were based. 
Information on the limiting of deceleration rates by modulation, for 
lifting vehicles, has not been available. 
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The purpose of the present report is to answer the following 
questions: Given a manned vehicle with a maximum lift -drag ratio of 0.5 

entering the earth’s atmosphere at a given flight -path angle , how much 
can the deceleration and rate of deceleration he reduced by modulation, 
and what are the corresponding changes in convective and radiative heating 
and range? To answer these questions a trajectory analysis has been made 
on an IBM 704 computing machine to solve the two-dimensional equations of 
motion of reference 8 . The results obtained are for parabolic entry 
velocity, a nonrotating spherical earth, and an exponential atmosphere. 

The results are applicable to vehicles that can be operated in the angle - 
of -attack range from maximum lift coefficient to minimum drag coefficient. 


NOTATION 

A reference area for drag and lift, ft 2 
C-p drag coefficient, ■ 

d) o drag coefficient at zero lift 

C T lift coefficient, . — . 

L pV 2 A 

D drag force , lb 

g local gravitational acceleration, ft sec -2 

G deceleration in g units 

L lift force, lb 

m mass of vehicle, slugs 

q total heat absorbed at the stagnation point, Btu ft"” 2 

4 heating rate at the stagnation point, Btu ft“ 2 sec _:L 

dt 

r distance from the center of planet , ft 

r G radius of planet, 2. 0926x10 r ft for earth 

R radius of curvature of vehicle surface, ft 

s range , ft 

t time, sec 
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u tangential velocity component normal to a radius vector, ft sec " 1 

v radial velocity component, ft sec " 1 

V resultant velocity, ft sec " 1 

y altitude, ft 

a angle of attack of vehicle, deg 

0 atmosphere density decay parameter, 1/23,500 ft -1 for earth 

7 flight -path angle relative to the local horizontal, negative for 
descent, deg 

p gravitational constant, 1.4078xi0 ls ft 3 sec” 2 for earth 
p atmosphere dens ity , slugs ft -3 

p Q atmosphere density at planet surface, O.OO238 slug ft" 3 for earth 

P 0 mean value for exponential approximation to atmosphere de ns ity - 
altitude relation, 0.0027 slug ft" 3 for earth 

Subscripts 

c convective 

1 init ial 

max maximum 

p vacuum perigee 

r radiative 

ANAIXSIS 

Trajectory Equations 


A trajectory analysis has been made utilizing the solution of two- 
dimensional equations of motion for entries into an exponential atmosphere 
of a nonrotating spherical earth. The polar coordinate system with 
velocity components, aerodynamic forces, and flight -path angle is defined 
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in the sketch. The differential 
equations for the velocity in the 
radial and tangential directions 
are* respectively, (see ref. 8) 
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dt 
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and g is the local gravitational acceleration given by 


The constant p in equation ( 5 ) is the gravitational constant defined by 
Newton’s inverse square law of gravitational attraction. The differential 
equations employed for the altitude and range are, respectively, 


dy 

dt 


v 


( 6 ) 


ds _ 
dt 


u 


( 7 ) 


The differential equations used for the total laminar convective heat 
(after ref. 8) and the radiative heat (ref. 9) absorbed per unit area at 
the stagnation point are, respectively. 


d( 4c _ • _ 16,600 j~p~ / V 'f 

at “ qc ' Jfb te j 

dq r . f ( V * V ) 

— = q r = RIO ; 


( 8 ) 

( 9 ) 


where R is the radius of curvature of the vehicle surface at the stag- 
nation point, p is the local atmosphere density given by 


( 10 ) 
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u tangential velocity component normal to a radius vector, ft sec" 1 

v radial velocity component, ft sec" 1 

V resultant velocity, ft sec" 1 

y altitude, ft 

a angle of attack of vehicle, deg 

3 atmosphere density decay parameter, 1/23,500 ft -1 for earth 

7 flight -path angle relative to the local horizontal, negative for 

descent, deg 

P gravitational constant, 1 . 4078 xi 0 16 ft 3 sec -2 for earth 
p atmosphere density* slugs ft -3 

p Q atmosphere density at planet surface* O.OO238 slug ft" 3 for earth 

p Q mean value for exponential approximation to atmosphere density- 
altitude relation* 0.0027 slug ft" 3 for earth 

Subscripts 

c convective 

i init ial 

max maximum 

p vacuum perigee 

r radiative 

ANALYSIS 

Trajectory Equations 


A trajectory analysis has been made utilizing the solution of two- 
dimensional equations of motion for entries into an exponential atmosphere 
of a nonrotating spherical earth. The polar coordinate system with 
velocity components* aerodynamic forces* and flight -path angle is defined 
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in the sketch. The differential 
equations for the velocity in the 
radial and tangential directions 
are , respectively, (see ref. 8 ) 
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and g is the local gravitational acceleration given by 


The constant (j. in equation (5) is the gravitational constant defined by 
Newton’s inverse square law of gravitational attraction. The differential 
equations employed for the altitude and range are , respectively. 
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The differential equations used for the total laminar convective heat 
(after ref* 8) and the radiative heat (ref. 9 ) absorbed per unit area at 
the stagnation point are, respectively. 
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where R is the radius of curvature of the vehicle surface at the stag- 
nation point, p is the local atmosphere density given by 


P » P 0 e-Py 
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and V is the resultant velocity given by 

V = »/u 2 + v 2 (ll) 

The radiative heating rate per unit area at the stagnation point is 
obtained from an interpolation of a table for the logarithm of equa- 
tion ( 9 ) f that is f 


l°gio(~) = f(y,v) ( 12 ) 

Values of (q r / R ) as a function of altitude and velocity for air in equi- 
librium were obtained from reference 9 . 


The six equations (l), (2), ( 6 ), ( 7 ), ( 8 ), and ( 9 ) were programed 
for simultaneous solution on an IBM 704 computing machine for entries 
during which the vehicle aerodynamic coefficients were constant and for 
entries during which these coefficients were varied (modulated) in order 
to satisfy specified constraints of a maximum deceleration rate and/or a 

maximum deceleration. The modulation started at Ct and ended at On 

■“max -“o 

because this type of modulation provides the lowest possible maximum 
deceleration (see, e.g., ref. 6). In the analysis the resultant decelera- 
tion G is given by 


G 


Pq 

2p(m/A) 


r 2 e -(3y V 2 C D 



( 13 ) 


For entries modulated to maintain a specified maximum deceleration rate, 
the deceleration time history during the modulation period is linear, that 
is , 


— — = constant ( l4 ) 

dt v ' 

and for entries modulated to maintain a specified maximum deceleration, 
the deceleration is constant during the modulation period. 


Vehicle Characteristics 


Modulation is accomplished in the present investigation by varying 
the vehicle drag coefficient and lift-drag ratio during entry. The drag 
coefficient and lift-drag ratio are calculated on the assumption (as in 
ref. 4 ) that the vehicle has a variation of lift and drag similar to that 
for a flat plate in Newtonian flow given by 


6 


c d - c d 0 + (yw 
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where 


h 



The particular drag polar with a maximum lift -drag ratio of 0.5 used in 
the present investigation is given in figure 1. The results , however, 
can he shown to apply to a family of vehicles having (L/D) max =0.5 with 
various values of *lD max and Cj) q provided Cp ma / Cp Q remains constant 

and the initial m/A is adjusted to give initial values of m/c^A equal 

to those used in the present investigation. A value of m/A = 3 was used 
in the present calculations . 


The results of the present analysis can he applied to vehicles of 
arbitrary weight and size hy employing the results of reference 8. Accord- 
ing to these results, many of the trajectory parameters are essentially 
independent of m/CpA. In this category are deceleration, deceleration 
rate, velocity, flight -path angle, and range. Trajectory parameters that 
depend on m/C-pA are altitude and convective and radiative heating where 
the relationships for altitude and convective heating can he shown to he 


and 


J, 


*1 



~ (m/c D A) x ~ 

_(m/c D A) 2 _ 



/(rn/OpAR)" 

J(m/0jyAR) x 


(18) 


(19) 


The subscript 1 refers to values of the present report and the sub- 
script 2 corresponds to other values of m/CpA or m/CpAR. Equations (9) 
and (l8) and reference 9 can be used to calculate radiative heating if 
the velocity -time and velocity-altitude relationships are given. 


A 
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Corridor Depth 


In. general, corridor depth is defined as a difference in perigee 
altitudes of vacuum trajectories corresponding to a difference in two 
flight-path angles at a given initial altitude and velocity. The flight- 
path angle associated with the higher of the two perigee altitudes is 
usually referred to as the angle for an overshoot "boundary, and the angle 
associated with the lower of the two perigee altitudes is referred to as 
the angle for an .undershoot "boundary (see, e.g., ref. ■7). The relation- 
ship "between perigee altitude and flight -path angle for a vehicle at 
parabolic velocity can be shown from vacuum trajectory relationships to be 

y p = y ± cos z y 1 - r 0 (l - cos 2 7i) (20) 

Perigee altitude and corridor depth are presented as functions of flight- 
path angle in figure 2 for an initial altitude of 400,000 feet. A corridor 
depth of zero was chosen so that it corresponded to 7 = -5.O3 0 . This is 
defined as the overshoot boundary in the present report since the vehicle 
entering at its highest negative lift-drag ratio (-0.5) skips out of the 
atmosphere for entries with angles of descent less than 5.03 0 . 


RESULTS AND DISCUSSION 


In order to determine the effects of modulation of lift and drag for 
a vehicle entering the earth's atmosphere at parabolic velocity, unmodu- 
lated and modulated trajectories have been computed. Unmodulated entries 
are defined in this report as entries during which the lift -drag ratios 
are held constant until the flight path is essentially horizontal (7=0, 
near maximum deceleration). To prevent skipping out of the atmosphere 
for the unmodulated entries with l/d = 0.5* step changes to l/D = -0.5 
then back to L/D = 0.5 were programed to occur at points along the flight 
path so that decelerations greater than the first peak decelerations were 
avoided. To prevent excessive decelerations for the unmodulated overshoot 
entry with L/D = -0.5, a step change to l/d = 0, a = 90°, was made at 
1 0* For the modulated trajectories, the initial lift-drag ratio was 

0.42, the value for maximum lift coefficient. Two types of modulation 
were employed; one held the maximum deceleration to a specified value, 
the other held the maximum deceleration rate to a specified value, with 
a further stipulation that a maximum deceleration of lOg was not to be 
exceeded. 

In figure 3 "typical modulated ’trajectories are compared with an 
unmodulated trajectory having identical initial conditions. Results 
from these and similar trajectories have been compiled to provide a 
summary of the effect of modulation on maximum deceleration and 
deceleration rate, total and peak stagnation point heating, and range. 
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Maximum Deceleration 


Figure 4 shows the maximum decelerations for unmodulated and modulated 
entries as functions of corridor depth and initial flight-path angle. For 
the 'unmodulated entries the maximum decelerations increased with increased 
corridor depth and decreased with increased lift -drag ratio as has teen 
found previously (e.g., refs. K, 5, and 8). For fixed corridor depth, the 
use of modulation provided reductions of maximum deceleration to values 
that were approximately 60 percent of the lowest unmodulated values 
(l/d = 0.5)* A few trajectories were also computed for initial velocities 
from 15,000 to 50,000 feet per second and the same reduction in maximum 
deceleration was found . For entries at escape velocity with a maximum 
deceleration limit of lOg, the above reduction in maximum deceleration is 
equivalent to lowering the undershoot boundary 2b miles . Thus , if the 
vehicle were operated at l/d = -0.5 on the overshoot boundary and modula- 
tion were used to enter at the steepest possible angle and not exceed lOg, 
the resulting corridor depth would be 62 statute miles . 

It is of interest to note that if modulation is attempted to reduce 
the maximum deceleration to values below those of the lowest curve of fig- 
ure b } a maneuver might be required after reaching a = 0° to prevent maxi- 
mum decelerations greater than that experienced if modulation had not been 
used. This is illustrated in figure 5 which shows the maximum deceleration 
obtained by modulation as a function of the maximum deceleration desired 
by modulation for an initial flight-path angle of -8.l4° (corresponding to 
a corridor depth of 50 miles). The maximum deceleration for an unmo dulated 
entry -under these conditions with l/d = 0.5 was about 13 . 2g. It can be 
seen from figure 5 that if modulation is attempted to hold the maximum 
deceleration to less than about 7g, the maximum deceleration will exceed 
the maximum deceleration for the unmodulated entry. This occurs because 
the rate of decrease of angle of attack required is so great that the 
flight -path angle is still relatively steep when the vehicle is operating 
at 0° angle of attack. Decelerations higher than those for the unm odulated 
case can be avoided, however, by an increase in angle of attack after 
a = 0° is reached. 


Maximum Deceleration Rate 


Figure 6 shows the maximum rates of increase of deceleration for 
unmodulated and modulated entries as functions of corridor depth and entry 
angle. For corridor depths up to about Mi miles, modulation was used to 
reduce the deceleration rates as much as possible without regard to the 
maximum decelerations. For corridor depths from 44 to 62 miles, modulation 
was used to reduce the deceleration rates as much as possible and still 
limit the maximum deceleration to lOg. For a given corridor depth, the 
deceleration rates were generally reduced to less than half the values for 
the unmodulated entries with l/d = (L/D) max =0-5- It can also be seen 
from figure 6 that an undershoot boundary that is limited by a maximum 
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deceleration of lOg and a maximum deceleration rate of 0.5g per second 
(38 mile corridor depth) can he lowered 20 miles (58 mile corridor depth) 
by modulation. This is only 4 miles (62 mile corridor depth) less than 
the undershoot boundary limited by a maximum deceleration rate of over 
lg per second. 


Stagnation-Point Heating 


Maximum convective h eati ng rate .- The maximum laminar convective 
heating rates for the unmodulated and modulated entries are shown in 
figure 7 as functions of corridor depth and entry angle. The heating 
rates generally increased with increasing corridor depth and entry angle. 
For a given corridor depth, the heating rates were approximately the same 
for . the unmodulated entries with L/D = 0.5 and the entries modulated to 
limit the deceleration; modulation to limit the deceleration rate generally 
resulted in heating rates about 20 percent higher than those for the unmod- 
ulated L/D =.0.5 entries. Modulation to lower the lOg-limited undershoot 
boundary 24 miles (corridor depth increase from 38 to 62 miles) resulted 
in an l 8 -percent increase in peak convective heating rate. 

Maximum ra diative heating rate .- The maximum radiative heating rates 
for the unmodulated and modulated entries are presented in figure 8 as 
functions of corridor depth and entry angle. The radiative heating rates 
generally increased with increasing corridor depth and entry angle. How- 
ever, for a given corridor depth, the peak radiative heating rates for the 
unmodulated entries with l/d = 0.5 were considerably higher than the rates 
for the entries modulated to limit the maximum deceleration. Modulation 
to lower the lOg -limited -undershoot boundary 2b miles resulted in about a 
10 -percent increase in peak radiative heating rate. 

Total convective heating .- The effect of modulation on the total 
convective heat absorbed per unit area at the stagnation point is shown 
in figure 9* To obtain the total heat for the unmodulated entries with 
L/D = 0.5, the vehicle was maneuvered as described earlier. Trajectories 
were also calculated for other maneuvers after peak deceleration to keep 
the vehicle from skipping out of the atmosphere ; for these maneuvers the 
total heat was about the same or greater than that shown. It can be seen 
from figure 9 that modulation to lower the lOg undershoot boundary 2b miles 
resulted in a decrease in total convective heating of 22 percent. 


Total ra diative heating .- It can be seen from figure 10 that modula- 
tion to lower the lOg undershoot boundary 2b miles resulted in about a 
10-percent increase in total radiative heating. Modulation to limit the 
deceleration rate for initial flight-path angles such that G max < 10 , 
however, resulted in much larger increases in total radiative heating. 
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Range 


Xt is evident from figure 11 that the range increased very rapidly 
for the unmodulated and modulated entries as the corridor depth or flight - 
path angle was decreased to such values that the vehicle was about to skip 
out of the atmosphere . Modulation to lower the lOg undershoot "boundary 
24 miles resulted in a decrease in range of about 60 percent. 


CONCLUDING REMARKS 


Unmodulated and modulated trajectories have been computed for a 
vehicle having a maximum lift-drag ratio of 0.5 entering the earth’s 
atmosphere at parabolic velocity. The results indicate that modulation 
can reduce the maximum deceleration to 60 percent of the unmodulated value 
or reduce the maximum deceleration rate to less than 50 percent of the 
unmodulated rate. These results are equivalent, for a maximum decelera- 
tion of lOg, to lowering the undershoot boundary 2 k miles with a resulting 
decrease in total convective heating to the stagnation point of 22 percent. 
However, the maximum convective heating rate was increased 18 percent.; the 
maximum radiative heating rate and total radiative heating were each 
increased about 10 percent. 
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7i = ^00,000 ft, m/A = 3 . 


obtained by modi 



G mQX desired by modulation 

Figure 5 .- Maximum deceleration obtained versus maximum deceleration desired: V- = ^6 335 ft/ 

y± = ^00,000 ft, 7j_ = -8.l4° (50 mile corridor), m/A = 3. 1 

» -p *. » . 


-f^O\Ul > 




CJ 


400 


-5 



10 



Unmodulated 

Modulated for minimum G 


m 


Modulated to limit ( 


ond G moxi'° 


* mo: 


20 30 40 50 60 

Corridor depth, statute miles 


-7 -8 

Initial flight -path angle, y , deg 


Figure 7«- Maximum laminar convective heating rates for unmodulated and modulated entries 
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